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Abstract 

We present a general form for the solution of an expanding general- 
relativistic Friedmann universe that encounters a singularity at finite fu- 
ture time. The singularity occurs in the material pressure and acceleration 
whilst the scale factor, expansion rate and material density remain finite 
and the strong energy condition holds. We also show that the same phe- 
nomenon occurs, but under different conditions, for Friedmann universes 
in gravity theories arising from the variation of an action that is an arbi- 
trary analytic function of the scalar curvature. 

There has been considerable interest in the possibility of finite and infinite 
future-time cosmological singularities in general relativity under the headings of 
'big-rip' singularities and associated 'phantom matter' evolution J,^. These 
studies have also extended to several quantum and brane-world generalisations 
13 , ^inii in]; miHI' a demarcation of these and related finite-time singularities 
has been given in [5]. 

We consider the most general form of sudden, or finite-time, singularity of 
the sort discussed in |1L)| that can arise in an Friedmann expanding universe. 
In jn] we showed that singularities of this sort can arise from a divergence of 
the pressure, p, at finite time despite the scale factor, a(t), the material density, 
p, and the Hubble expansion rate, H = a/a remaining finite. The pressure 
singularity is accompanied by a divergence in the acceleration of the universe, 
a, at finite time. Remarkably, these singularities occur without violating the 
strong energy conditions p > and p + > 0. They can even prevent a 
closed Friedmann universe that obeys these energy conditions from attaining an 
expansion maximum In order to prevent the occurrence of these sudden 

singularities it is necessary to bound the pressure. A sufficient condition is to 
require dp/ dp to be continuous or p/p to be finite jlOj . 

Consider the Friedmann universe with curvature parameter fc, and 8nG = 
c = 1, the Einstein equations reduce to 
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3i?2 = p-A, (1) 

= p + iH{p + p), (2) 

1 - -(^)- (^) 

In we constructed an explicit example by seeking, over the time interval 
< i < tg, a solution for the scale factor a{t) of the form 

a(t) = l + Bt'i + C(ts-tY\ (4) 

where i? > 0, g > 0, C and n> Q are free constants to be determined. If we fix 
the zero of time by requiring a(0) = 0, so Ci" = —1, we have 

«w = (f)M«s-i) + i-(i-fr, (5) 

where Og = a{ts)- Hence, as t ^ from below, we have 

a^^li'l-DBf-^-^^j^^-^ (6) 

whenever 1 < n < 2 and < q < 1; the solution exists on the interval < t < t^. 
Hence, as t —^ ts we have a — > Og, iJ — > and p ^ > where as, -ffs, and 
Ps are all finite but ps —^ oo. By contrast, as t ^ there is an initial all- 
encompassing strong-curvature singularity, with iJ ^ oo, p — > oo and p — > oo. 
From © and we see that p and p + 3p remain positive throughout the 
evolution but, because p is finite asymptotically, the dominant-energy condition, 
\p\ < p, must always be violated as i — > tsj 10 . A solution with similar 
properties which expands from a de Sitter past state also exists, with 

a{t) = fls - 1 + exp{A(ts - t)} - (1 - ■^)" (7) 

ts 

with A > constant and 1 < n < 2. 

The particular solution ^ expands with 



as i ^ and so with choices q = 1/3,1/2, 2/3, or g = 1 it will resemble a 
massless scalar-, radiation-, dust-, or negative-curvature-dominated Fricdmann 
universe respectively at early times. At late times, as the 'big-rip' singularity is 
approached, the expansion approaches a constant with 

a{t) ^ as + q{l - as){l - ^) (9) 

as t ^ to ■ 
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Consider now the most general form of a solution of this type. In the neigh- 
bourhood of the singularity at t = ts the general formal solution of the Fried- 
mann equations has the form 



a{t) = {t, - t)" 



Ni 



j=0 k=0 



(10) 



where [...] is a convergent double Psi-series which tends to zero as i — > tg; Ujk 
are constants, Nj < j are positive integers and Q G (this is derived as an 
application of the theorem of Goriely and Hyde proved in ref. p.2J). 
Hence, we have a generalisation of the solution © of the form 



a{t) = i-y (as - + its -tr 



Ni 



j=0 k=0 



(11) 



As t ^ ts if we choose 1 < rt < 2, < (7 < 1, we have a 



and 



d -nin - l)A{ts - t)""^ log''(ts - t) + higher order -00, (12) 

where A > is constant, and so p + 3p > and p > \p\ as the future singularity 
is approached with p ^ 00 at ts- As i the solution evolves as a w cq + cif 
and approaches the behaviour of the radiation-dominated Friedmann universe 
for q = l/2. Ast^oo the solution approaches Ug in accord with eq.©. We 
see that the solution © is the special case with all ajk = except for aoo 7^ 0. 

By choosing n to lie in the intervals {N,N + 1) for > 2 where N e Z+, 
we can create a singularity in which 



but 



dt^+^ 



0, forr<Ne Z+. 



(13) 



(14) 



In particular, this allows for singularities of the pressure which are accompa- 
nied by divergences of higher time-derivatives of a(t) in Friedmann solutions of 
higher-order gravity theories • Consider a generalisation of general relativ- 
ity to a theory of gravity derived from the variation of a lagrangian that is an 
arbitrary polynomial function /(i?) of the scalar curvature, R. The Friedmann 
universe is now prescribed by the solution of any two of equations ^ and the 
generalisations of Q and which are given respectively by ^21 



Rf - 2/ + 3f"{R + S-R) + if'R^ = Sp-p 
a 



(15) 



3 



and 



a , 1 ,, ■ a , ^ 

3-f' + -f-3f'R-+p = 16 
a Z a 

where /' = df / dR etc, and the scalar curvature of the Friedmann metric is given 

by 

^ „ f d k\ ,^ „. 

i?- -6 - + - + - . 17 

The equations of general relativistic cosmology in the absence of a cosmological 
constant arise from H15|l - H17|) as the special case f{R) = R. Various auxiliary 
conditions on the form of f{R) are required to ensure the existence of Friedmann 
solutions, see ref. 13 for a discussion. We ignore them as they do not affect 
our results here. 

In order to obtain a finite-time singularity of the sort described by |(SJ) or 
pi|l we need to accommodate a singularity of the pressure, p, in |(2Jl in equa- 
tions (I15|) - (I16|I . We require that all the terms in the generalised Friedmann 
first integral l|l()|l remain finite as t — > and a — > Os. Thus we require that 
a, f{R), f'{R)i f"{R), R, R and p are all finite in this limit. 

Now consider the evolution equation (|15|l a.s t ^ tg. Since R,f,f',f" and 
/"' will be finite at ts we must have 

f'R^p ^ -oo as t^ts. (18) 
Hence, there must be a divergence of 

R -6— -oo ast^ts (19) 
a 

but with R{ts) and R{ts) remaining finite. This occurs so long as 3 < n < 4, 
since for a solution of the form © 

hm ^ = _n(n-l)..(n-r+l) _ ^ 

t^t, dr % ts 

and the same condition holds for the general solution Hll|l . In the general 
relativity case of f{R) — R, the pressure singularity accompanies one in d and 
requires 1 < n < 2. 

This generalisation of the conditions for the existence of sudden singularities 
to the Friedmann universes of higher-order lagrangian theories of gravity shows 
that the higher-order curvature corrections that these extensions introduce are 
able to moderate but not remove the finite-time singularity. They merely change 
the evolutionary behaviours of the expansion scale factor which permit their 
occurrence. As in general relativity, an energy condition that bounds p/p 
suffices to exclude them. 



4 



References 



[1] Starobinsky A A 2000 Grav. Cosmol. 6 157; Caldwell, R R 2002 Phys. Lett. 
B 545 23; Shulz A and White M J 2001 Phys. Rev. D 64 043514; Hao J 
and Li, X 2003 Phys. Rev. D 67 107303; Gibbons G W hep-th/0302199; 
Nojiri S and Odintsov S D 2003 Phys. Lett. B 562 147 and B 571 1; Singh 
P, Sami M and Dadhich N 2003 Phys. Rev. D 68 023522; Hao J and Li, 
X 2003 Phys. Rev. D 68 04351 and 083514; D§,browski M, Stachowiak T 
and Sydlowski M 2003 Phys. Rev. D 68 067301; Elizalde P and Quiroga J 
2004 Mod Phys Lett A 19 29; Gonzalez-Diaz P F 2004 Phys Lett B 586 
1; Feinstein A and Jhingan S 2004 Mod Phys Lett A 19 457; Chimento L 
P and Lazkoz R 2003 Phys Rev Lett 91 211301; Elizalde E, Nojiri S, and 
Odintsov S D 2004 Phys. Rev. D 70 043539; Chimento L P and Lazkoz R 
2004 Mod Phys Lett A 69 123512. 

[2] Gonzalez-Dfaz P F 2003 Phys. Rev. 68 021303(R) 

[3] Calcagni G 2004 Phys Rev D 69 103508 

[4] Gorini V, Kamenshchik A, Moschella U and Pasquier V 2004 Phys Rev D 
69 123512 

[5] Kofinas G, Maartens R and Papantonopoulos E 2003 JHEP 0310 066 

[6] Sahni V and Shtanov Y 2002 Class. Quantum Grav. 19 LlOl 

[7] Nojiri S and Odintsov S D 2004 Phys. Lett. B 595 1 

[8] Nojiri S and Odintsov S D Preprint hep-th/0408170] 

[9] Cotsakis S and Klaoudatou I Preprint gr-qc/040902^ 

[10] Barrow J D 2004 Glass. Quantum Grav. 21 L79 

[11] Barrow J D, Galloway G J and Tipler F J 1986 Mon. Not. R. astron. Soc. 
223 835 

[12] Goriely A and Hyde C 2000 J. Dijf. Eqns. 161 422 

[13] Barrow J D and OttewiU A C 1983 J. Phys. A 16 2757 



5 



